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Abstract. Using a new radio sample, 6C* designed to find radio galaxies at z > 4 
along with the complete 3CRR and 6CE sample we extend the radio galaxy K — z 
relation to z ~ 4.5. The 6C* K — z data significantly improve delineation of the K — z 
relation for radio galaxies at high redshift (z > 2). In a spatially flat universe with 
a cosmological constant (17 m = 0.3 and 11 a = 0.7), the most luminous radio sources 
appear to be associated with galaxies with a luminosity distribution with a high mean 
(« 5 L*), and a low dispersion (a ~ 0.5 mag) which formed their stars at epochs 
corresponding to z J 2.5. 


1 The advantages of using radio galaxies 

Radio galaxies provide the most direct method of investigating the host galaxies 
of quasars if orientation based unified schemes are correct. The nuclear light 
which dominates the optical/near-infrared emission in quasars is obscured by 
the dusty torus in radio galaxies, therefore difficult psf modelling and subtrac¬ 
tion are not required to determine the properties of the underlying host galaxy. 
Unfortunately compiling samples of radio loud AGN is a long process, because 
of the radio selection there is no intrinsic optical magnitude limitation, making 
follow-up observations extremely time consuming, especially when dealing with 
the faintest of these objects. However, low-frequency selected radio samples do 
now exist with the completion of 3CRR (Laing, Riley & Longair 1983) along 
with 6CE (Eales et al. 1997; Rawlings et al. 2001) and the filtered 6C* sample 
(Blundell et al. 1998; Jarvis et al. 2001a; 2001b). We can now use these ra¬ 
dio samples to investigate the underlying stellar populations through the radio 
galaxy K — z Hubble diagram. 

2 Previous radio samples and the K — z Hubble diagram 

There has been much interest in the K — z relation for radio galaxies in the past 
decade. Dunlop & Peacock (1993) using radio galaxies from the 3CRR sample 
along with fainter radio sources from the Parkes selected regions demonstrated 
that there exists a correlation between radio luminosity and the A"—band emis¬ 
sion. Whether this is due to a radio luminosity dependent contribution from a 
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non-stellar source or because the galaxies hosting the most powerful radio sources 
are indeed more massive galaxies has yet to be resolved. Eales et al. (1997) con¬ 
firmed this result and also found that the dispersion in the A'—band magnitude 
from the fitted straight line increases with redshift. This result, along with the 
departure to brighter magnitudes of the sources at high redshift led Eales et al. 
to conclude that we are beginning to probe the epoch of formation of these mas¬ 
sive galaxies. Using the highest redshift radio galaxies from ultra-steep samples 
of radio sources van Breugel et al. (1998) found that the near infrared colours of 
radio galaxies at z > 3 are very blue, consistent with young stellar populations. 
They also suggest that the size of the radio structure is comparable with the size 
of the near infrared region, and the alignment of this region with the radio struc¬ 
ture is also more pronounced at z > 3. Lacy et al. (2000) using the 7C-III sample 
found evidence that the hosts of radio galaxies become more luminous with red¬ 
shift and are consistent with a passively evolving population which formed at 
high redshift (z > 3). Thus, all of this work points to a radio galaxy population 
which formed at high redshift and has undergone simple passive evolution since. 
However, all of these studies were made with only a few high-redsliift (z > 2) 
sources. With the 6 C* sample we are now able to probe this high-redshift regime 
with increased numbers from samples with well-defined selection criteria. 


3 The 6C* filtered sample 

The 6 C* sample is a low-frequency radio sample (0.96Jy < S 151 < 2.00 Jy) which 
was originally designed to find radio sources at z > 4 using filtering criteria 
based on the radio properties of steep spectral index and small angular size. The 
discovery of 6C*0140+326 at z = 4.41 (Rawlings et al. 1996) and 60*0032+412 
at z = 3.66 (Jarvis et al. 2001a) from a sample of just 29 objects showed that 
this filtering was indeed effective in finding high-redshift objects. Indeed, the 
median redshift of the 6 C* sample is z ~ 1.9 whereas for complete samples at 
similar flux-density levels the median redshift is z ~ 1.1. We can now use this 
sample to push the radio galaxy K — z diagram to high redshift (z > 2) where 
it has not yet been probed with any significant number of sources. Fig. |l| shows 
the radio luminosity-redshift plane for the three samples used in this analysis. 

4 Emission-line contamination 

The most-luminous sources at high redshift may be contaminated by the bright 
optical emission lines redshifted into the infrared. This is particularly true for 
sources in radio flux-density limited samples. The high-redshift sources in these 
samples are inevitably some of the most luminous, and we also know there is 
a strong correlation between low-frequency radio luminosity and emission-line 
strength (e.g. Rawlings & Saunders 1991; Willott et al. 1999; Jarvis et al. 2001a) 
which will increase the contribution to the measured AT—band magnitudes from 
the emission lines in the most radio luminous sources. 
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Fig. 1. Rest-frame 151 MHz luminosity (T 151 ) versus redshift 2 plane for the 3CRR 
(circles), 6 CE (triangles) and 6 C* (stars) samples. The rest-frame 151 MHz luminosity 
Z /151 has been calculated according to a polynomial fit to the radio spectrum (relevant 
radio data from Blundell et al. 1998). The curved lines show the lower flux-density 
limit for the 3CRR sample (dotted line; Laing et al. 1983) and the 7CRS (solid line; 
Blundell et al. in prep; Willott et al. in prep). The dashed lines correspond to the 
limits for the 6 CE sample (Rawlings et al. 2001) and the shaded region shows the 6 C* 
flux-density limits (all assuming a low-frequency radio spectral index of 0.5). Note that 
the area between the 3CRR sources and 6 CE sources contains no sources, this is the 
area which corresponds to the absence of a flux-density limited sample between the 
6 CE (S 151 < 3.93 Jy) and 3CRR (Sirs > 10.9 Jy) samples. The reason why some of the 
sources lie very close to or below the flux-density limit of the samples represented by 
the curved lines is because the spectral indices lie very close to or below the assumed 
spectral index of the curves of a = 0.5. 


To subtract this contribution we use the correlation between [Oil] emission¬ 
line luminosity T[oii] and the low-frequency radio luminosity L 151 from Willott 
(2000), where T[oii] cx Z'i 5 0 l 0±0 ' 04 . Then by using the emission-line flux ratios for 
radio galaxies (e.g. McCarthy 1993) we are able to determine the contribution to 
the A'—band magnitude from all of the other emission lines. This is illustrated 
in Fig. H where the emission-line contamination to the AT—band flux is shown 
for various radio flux-density limits and a range of redshifts. 


5 The K — z relation 

In Fig. 3 we plot the K — z diagram for all of the sources in our dataset. We 
also show four synthetic galaxy evolution models from the ‘Galaxy Isochrone 
Synthesis Spectral Evolution Library’ (GISSEL) of Bruzual & Chariot (1993) 
and a curve representing a galaxy which undergoes no-evolution. The GISSEL 
files that we have used are ones in which there is an instantaneous burst of 
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Fig. 2. Emission line contribution to the A'—band magnitudes for various radio flux- 
densities assuming the power-law relation of L[ 0 n| oc L\'Si- 


star formation and one in which the burst of star formation lasts 1 Gyr, with 
a Salpeter IMF with a lower mass cut-off of 0.1 M 0 and an upper mass cut-off 
of 125 M©. We use two different assumptions about the star formation history, 
one in which the burst of star formation begins at 2 = 5 and one in which the 
burst occurs at z = 10. 

The no-evolution curve was constructed by taking the spectral energy dis¬ 
tribution template from the GISSEL library that was found to fit the observed 
spectral energy distribution of a radio galaxy at 2 = 0, and which also repro¬ 
duced the near-infrared colours. All of the curves are normalised to pass through 
the low-redsliift (z < 0.3) points. 

With our data on the 6C* sample in addition to the 6CE and 7C-III samples 
we find that in a low-density universe (Qm = 0.3 and = 0.7) the data 
are predominantly brighter than the no-evolution curve and are consistent with 
a passively evolving stellar population with a high-formation redshift. If this 
passively evolving scenario is correct then hierarchical growth at z < 2.5 is not 
a required ingredient. However, this brightening may not just be due to passive 
evolution of the stellar population. Non-stellar contributions from the central 
AGN may also contribute a higher fraction of light at these redshifts. All of the 
studies to measure the non-stellar contribution to the A'—band flux conducted 
to date (e.g. Leyshon & Eales 1998; Simpson, Rawlings & Lacy 1999), have 
concentrated on the most radio luminous 3CR sources at z ~ 1, and may have 
little bearing on the properties of the high-redshift sources considered here. If 
it turns out that the high-redshift 6C sources have > non-stellar contamination 
to those of the 2 ~ 1 3C sources (which have the same radio luminosity) then 
hierarchical build up may be necessary. Note that K— band observations of the 
high-redshift 6C sources will be at shorter wavelengths than those of the 3C 


sources. 
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Fig. 3. The K — z Hubble diagram for radio galaxies from the 3CRR (circles), 6CE 
(triangles), 6C* (stars) and 7C-III (squares) samples. A' 63.9 denotes the A'—band mag¬ 
nitude within a comoving metric aperture of 63.9 kpc (c.f. Eales et al. 1997; Jarvis et 
al. 2001b). The two vertical lines show the redshift above which the alignment effect 
begins to be seen (z > 0.6) and the higher redshift at which we chose to split the data 
(z = 1.8). The solid curved line is the predicted curve for galaxies which do not evolve 
(as described in the text). The dashed lines are models for a star-burst lasting 1 Gyr 
starting at z = 5 (lower) and z = 10 (upper). The two dotted curves represent the 
models of an instantaneous (0.1 Gyr) burst beginning at z = 5 (lower) and z = 10 
(upper). 


However, separate arguments lead us to conclude that the dominant factor 
is passive evolution of a stellar population which formed at z > 2.5. First, recent 
sub-nun observations with SCUBA have shown that the dust masses in radio 
galaxies are larger at z ~ 3 than in galaxies with similar radio luminosities at 
lower redshift (Archibald et al. 2001). This implies that the majority of star- 
formation activity in these galaxies is occurring at high redshift. Second, the 
discovery of six extremely red objects at 1 < z < 2 in the 7C Redshift Survey 
(Willott, Rawlings & Blundell 2001) with inferred ages of a few Gyrs, implies 
that these objects formed the bulk of their stellar population at z ~ 5. Third, de¬ 
tailed modelling of the optical spectrum of the weak radio source LBDS 53W091 
at z = 1.552 has shown that this object is most plausibly an old elliptical, with 
an inferred age of >3.5 Gyr (Dunlop et al. 1996). The further discovery of LBDS 
53W069 at z = 1.43, with an inferred age of ~ 4 Gyr (Dunlop 1999) suggests 
that there exists a population of evolved, radio weak ellipticals which formed 
at z>5. Therefore, the new data on the 6C* sample presented in here is consis- 
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tent with the results from various other observational studies of radio galaxies 
in which these radio-luminous systems formed most of their stars at epochs cor¬ 
responding to very high redshifts (z >2.5), and have undergone simple passive 
stellar evolution since. Willott et al. (2001) have pointed out that such galax¬ 
ies probably undergo at least two active phases: one, at epochs corresponding to 
z>5, when the black hole and stellar spheroid formed, and another, at e.g. z ~ 2, 
when powerful jet activity is triggered, or perhaps re-triggered, by an event such 
as an interaction or a merger. The small scatter in the K — z relation (Jarvis et 
al. 2001b) and sub-mm results (Archibald et al. 2001) suggest that the second 
active phase has little influence on the stellar mass of the final elliptical galaxy 
We have shown that the powerful radio galaxies in our samples are consistent 
with having passively evolving stellar populations. If we now compare the masses 
of these powerful radio galaxies to the derived value of M^ for nearby ellipti¬ 
cal galaxies [M£- = —24.3 for H 0 = 70kms _1 Mpc _1 (Kochanek et al. 2000)], 
we find, if passive evolution is accounted for, that the powerful radio galaxies 
considered in this here are consistent with being sa 5 L* throughout the redshift 
range 0 < z< 2.5. 
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